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Abstract 



In this thesis the interaction of gravitational waves (gws) with electromagnetic 
waves (emws) in a static magnetic background field is considered. This interaction 
is a consequence of the general relativistic Einstein field equations. These equations 
dictate the excitation of electromagnetic waves when a gravitational wave interferes 
with a static electromagnetic background field and a similar reverse process. 

Such interactions can become effective, close to very energetic sources such as 
colliding neutron star binaries (i.e. gamma ray bursts), and quacking supernova 
remnants (magnetars), or in the large scale magnetic fields of the early universe. 
The former two sources have strong, localized and the latter weak, but extended 
magnetic fields. This is important because the energy transfer efficiency of GW <^> 
EMW conversions appears to be quadratically proportional to the background field 
strength and the extension of the interaction region. 

All calculations in this thesis are done in a general relativistic, space+time, non- 
coordinate formalism. The reason for this is that in such tedrad systems, equations 
remain transparent and facilitate easy physical interpretation and connection to 
measurements. 

In this framework, the conversion efficiency of GWs to light and vice versa in a 
vacuum is considered, first in an estimate and then in a more elaborate and gen- 
eral, exact calculation. The GW=^EMW conversion is proposed as a possible indirect 
detection device for gamma ray bursts and magnetars. Also, the possibility is con- 
sidered of explaining the small fluctuations in the cosmic background radiation by 
the conversion of Gws in the early universe to emws superposed on the homogeneous 
background radiation. 

Next, to obtain more realistic results, the same process is examined in a thin 
plasma which leads to the same emws as generated in a vacuum. The importance 
of the presence of this plasma, though, is that it might damp the generated radio 
waves before they can travel over astromical distances, unless non-linear effects lead 
to higher frequencies of the EMWs. If radio waves with large enough frequencies 
are generated, gamma ray bursts and supernova? might be detectable with (space 
based) radio detectors such as the proposed Astronomical Low Frequency Array 
(alfa) as well as with GW detectors such as Laser Interferometer Gravitational 
wave Observatory (ligo) with a event rate of as many as a few per year in our local 
galaxy group and the Virgo cluster. 

In the last chapter, an entirely different interaction is proposed, in which the 
gravitational wave interacts with the plasma and generates fast magneto-acoustic 
plasma waves, thus dissipating its energy into the plasma. The plasma can then 
emmit the energy as electromagnetic radiation. Theoretical models for gamma ray 
burst could be improved a lot if even a small fraction of the GW energy could be 
converted into emws in this fashion. The reason for this is that the energy released 
in a neutron star binary merger is expected to be released mainly in Gws, whereas 
the obeserved energy flux is mostly electromagnetic. The dispersion relation derived 
for this interaction is the most interesting, new, result of this thesis work. 
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1 INTRODUCTION AND OVERVIEW 



1 Introduction and overview 

A consequence of Einstein's General Theory of Relativity is the existence of the inter- 
esting concept of gravitational waves. These waves originate from the most energetic 
events in our universe, such as colliding neutron star binaries, supernova explosions and 
gravitational collapses into black holes. They manifest themselves as ripples in space and 
time, that periodically stretch and compress all present matter and delay and accelerate 
the time signals from millisecond pulsars, travelling to us over astronomical distances. 

1.1 Direct gravitational wave detection 

The only evidence, until now, of the existence of gravitational waves is the neutron 
star binary whose orbital radius is slowly decreasing as a consequence of the radiated 
gravitational energy. The reason that it is so hard to detect gravitational waves (gws) 
is that Newton's constant of gravitation is extremely small (G = 6.7 • 10~ n Nm 2 kg -2 ). 
The dimensionless amplitude (A///) of typical gravitational waves reaching the earth is 
only of the order of 10~ 17 . In other words, a rod of one meter in length will oscillate 
with an amplitude of a millionth of the radius of a hydrogen atom. Still, several direct 
GW-detectors are being build at present, or have been proposed for the future, that 
hope to detect these very small vibrations. The best known of these are the Laser 
Interferometer Ground Observatory, ligo, and the Laser Interferometer Space Array, 
LISA. As the names suggest, the former is a terrestial observatory and the latter a 
space base one. Needless to say wha enormous engineering achievements these detectors 
require to detect the oscillations mentioned above. 

Another approach to detect GWs is to measure 'relative time displacements' (At/t), 
instead of position displacements. Time signals from millisecond pulsars can be mea- 
sured to such precision that time delays caused by the gravitational waves can in principle 
be measured with the same precision as position displacements as soon as enough stable 
pulsars can be observed for a sufficiently long time. Dr. A. Lommen et. al. (Berkeley, 
U.S.A.) plan to use this fact to determine, from GWs, whether the central object of our 
galaxy is a black hole binary. 

1.2 Watching gravitational waves 

The subject of the first part of this thesis is an entirely different way to indirectly observe 
gravitational waves. It appears that GWs can be converted to electromagnetic waves, viz 
light, that are much easier to detect. 

According to General Relativity, all forms of energy are equivalent, which is one way 
to formulate the Equivalence Principle. This is reflected by Einstein field equations, 
which are the equivalent of the non-relativistic Poisson equation for the gravitational 
potential. Just as the latter gives the gravitational potential as a function of the matter 
density, the former gives the space-time curvature due to the total energy-momentum 
density. Gravitational waves are the first order, harmonic solutions for these field equa- 
tions. 
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1 INTRODUCTION AND OVERVIEW 



1.3 Cosmic background radiation 



As a result of the Equivalence Principle, not only oscillating matter sources can pro- 
duce GWs, but also oscillating em energy densities. In linearized theory, this means that 
also the reverse process might occur: electromagnetic waves generated by gravitational 
waves. It will be shown in Sec. |3| that a gravitational wave passing through a constant 
electromagnetic background field is partially converted into light. In other words, one 
could indirectly 'watch' gravitational waves in the form of light. 

In Sec. ^| an exact calculation will be presented for a linearly polarized emw passing 
through an arbitrary EM background field which converts it to a GW which is then 
reconverted into an EMW. From this calculation it is apparent that such conversions are 
absolutely inefficient in any laboratory experiment. From an observational view point, 
however, they could provide a means to indirectly observe the GWs from more energetic, 
astrophysical phenomena such as supernovas or merging neutronstar binaries. 



1.3 Cosmic background radiation 

As a second phenomenon where GW to EMW conversions could play an important role, the 
observed fluctuations in the cosmic background radiation are investigated. According to 
the standard model big bang theories, 20.000 years after the Big Bang, the electron mist 
that until then obscured the universe, evaporated, and the universe became transparent 
to the ~ 2700 Kelvin thermal radiation. 

This radiation is observed as a very isotropic flux of photons, coming to us from all 
directions. What has not been explained properly are the small fluctuations that do 
appear in the radiation. If, however, a primordial magnetic field existed in the early 
universe, gravitational waves travelling through this field could have been converted 
into electromagnetic waves causing small fluctuations in the otherwise homogeneous 



radiation. An estimate of this effect is given in Sec. 8.2. 



1.4 Gamma ray bursts 

A final GW to emw conversion is discussed in the last part of this thesis. In a strongly 
magnetized plasma surrounding a source of GWs, these waves might generate plasma 
waves, which in turn produce light in the form of cyclotron, synchrotron etc. radiation. 
In this case, the interaction is caused by the fact that in a very strong magnetic field, 
the plasma electrons are frozen to the field lines. The GWs cause displacements of the 
electrons, which thus induce oscillations in the magnetic field. 

The relation to gamma ray bursts is that in these bursts, most of the energy is 
expected to be released in the form of gravitational waves, whereas the observed energy 
flux is electromagnetic. A way to solve this problem would be to dissipate a small 
fraction of the gravitational wave energy into the plasma, which can then radiate this 
energy as the observed electromagnetic waves. This is a compelling alternative, for 
instance, to theories suggesting a neutrino fueled fireball. 
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1 INTRODUCTION AND OVERVIEW 



1.5 Overview 

The next section, Sec. 0, is a brief introduction to some of the concepts and equations 



that are needed in later sections. In Sec. 2.1 gravitational waves are derived as weak 



field solutions of Einstein's field equations and the non-coordinate tedrad formalism, 



which will be used for most of the calculations, is explained in Sec. 2.2. 

In Sees. [3f-|4], the efficiency of GW to emw conversion in a vacuum is considered, first 
in an estimate following Gertsenstein Jl5| in Sec. || and then, in an exact calculation, for 
an arbitrary em background field and arbitrary GW components in Sec. [|. Finally, as 
an astrophysical application, an order of magnitude calculation is given for GW to EMW 
conversion close to a quacking supernova remnant neutron star, or magnetar, in Sec. [El. 
In Sec. H the vacuum is replaced by a thin plasma leading to some dispersion of the 



generated light (section 6.4) and finally, in Sec. [7], pressure gradients are added and the 



magnetohydrodynamic approximation is considered. From this the plasma waves are 
derived, needed to provide a GW dissipation mechanism for gamma ray bursts. 

The thesis ends with some conclusions in chapter || (and in appendix [C] second order 
GW-effects are taken into account resulting in the excitation of longitudinal Alfven-like 
waves) . 



1.6 Note on units 



Throughout this thesis, Gaussian CGS units are used in all equations. In most exact 
calculations the speed of light and Newton's constant are suppressed by choosing c = 
G = 1 but in numerical estimates, these constants will sometimes reappear to emphasize 
their importance. 

Furthermore, in the first two chapters greek indices will conventionally indicate four- 
vector components, fj, = 0, 1, 2, 3, with time as the zeroth component, and latin indices 
the spatial components i,j = 1,2,3. In the last chapter, however, the indices are used 
the other way around, following the particular convention used by the authors of the 
articles discussed in that chapter and |27]-|pl|). 

Finally, the signature of flat-space metric is chosen as: 



A1V 



( 


-1 








o > 







1 
















1 





\ 











1 ) 



(1) 



JOACHIM MOORTGAT 



WATCHING GRAVITATIONAL WAVES 



7 



2 GENERAL RELATIVITY 



2 General relativity 

2.1 Gravitational waves 

2.1.1 Introduction 

Gravitational waves, just like any other type of waves, are defined as propagating per- 
turbations of some flat background. Just as water waves are considered to be small 
ripples on an otherwise flat ocean, GWs are identified as small ripples rolling across 
spacetime. And just as one neglects the deviation from flatness of the ocean caused by 
the curved surface of the earth, the tidal forces caused by the gravitational attraction 
of the sun and the moon, Coriolis forces due to the earths rotation etc., so one ignores 
the large-scale curved structure of space-time caused by the matter distribution or for 
instance the presence of a primordial magnetic background field. The GWs originating 
from supernovae, explosions in the galactic center, or rotating binary stars can then 
indeed be described as small ripples on a flat background, as illustrated by the cartoon 
picture from LISA, Figure |l], showing the gravitational waves of a binary. 

As a result of all this, the (linearized) wave equations are easy to derive, as will 
be done in this section. But one has to keep in mind that these equations are defined 
only locally, and that they have no meaning globally, where the large scale space-time 
structure of the universe comes into play. 

2.1.2 Linearized theory in vacuum 

In the General Theory of Relativity, Poisson's equation of gravity is replaced by the 
equivalent, but covariant field equations, introduced by Einstein: 

PoiSSON (2) 
Einstein (3) 

(4) 

The conceptual meaning of these equations is that a localized density distribution curves 
the space around it and as a result of this, the geodesies for light or particles are deflected 
in the direction of the center of mass, thus reproducing the effect of a gravitation force. 

To derive the gravitational wave equations, the deviations from the Lorentz metric 
{Viiv) are assumed to be small. In other words, the weak-field solutions of Einstein's 
field equations (||) are obtained by assuming: 



G 



or 



R 



R 



fJ,U 



^g^R 



□$ = 47175 

T 



1 



g^u = V^ + h^u where |/t Mi ,| <C 1 (5) 
K X = ? ? AQ V and h = h^ = rf x h aX 

which is called the linearized theory of gravity. In this metric the Ricci tensor to first 



order in h^ v is (see for instance [ 18 1 ) : 



WATCHING GRAVITATIONAL WAVES JOACHIM MOORTGAT 



2.1 Gravitational waves 



2 GENERAL RELATIVITY 



- r' 3 - 




- - 






a ~ 2^,^ ~ h ff ,v& ~ h" v> up) 


1, 
2 




= — 2^^' 


a 
a 



where in the last line, the coordinate conditions are chosen such that (h/f — \f}y?h\p = 
(i.e. a divergence free and traceless solution). The field equations are now given (to first 
order) by: 



!, rr 1 (7 SvrG^ 

where (f)" is defined by: 

K = K-\^ h with <.- = ° ( 8 ) 

The general solutions of (Q) are, as mentioned in the previous section, just the solutions 
of the Poisson equation: 

4G f (T,/) iet dV 
t- ./ r — r' 



2.1.3 Plane wave solutions 

In a vacuum and the absence of EM fields, where T^ v = 0, the field equations (|7j) reduces 
to h" — \r]^h = with as simplest solutions plane waves: 

with (10) 

k a k a = NULL VECTOR 
A^ u k U = TRANSVERSE 

These solutions are transverse plane waves propagating with the speed of light (u = |k|) 
in the direction \c/u. 

The amplitude of this wave still seems to have six independent components (ten 
components because A^ u is symmetric (since T^ u is) of which four are fixed by the 
restriction A iU ,k v = 0). This is due to a gauge freedom under infinitesimal coordinate 
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2.1 Gravitational waves 




transformations x'^ — > + £^ satisfying □£ = 0. In this case £^ = —iC^expik a x a 
generates gauge transformations that can arbitrarily alter four of A^s components, 
therefore one has to choose not only a certain coordinate condition (as in (||)) but also 
a specific gauge. 

2.1.4 Transverse traceless gauge 

Defining the observer velocity by u u , one can perform a gauge transformation such 
that A^ u u u = which fixes three of the components, and another transformation to set 
A^n = wich fixes the last one. What remains are the two dynamical degrees of freedom 
of the gravitational field (the two polarizations). Summarizing: 

V"" = A llv k v = A^ = Q (11) 
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(12) 



One result of this gauge is that cp^ = h^ u which makes the analogy with Poisson's 
equation complete. 

2.1.5 Polarization 

Just like electromagnetic waves, all gravitational waves can be decomposed into two 
linearly polarized components, or in two circularly polarized components that give the 
geodesic deviation in a certain direction (0]). Linearly polarized, electromagnetic waves 
with polarization vectors e x or e y , travelling in the z-direction, cause a test particle to 
oscillates in the x- or y-direction respectively (with respect to an inertial frame). 
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2 GENERAL RELATIVITY 



For a gravitational wave, from (|l2|), the unit linear-polarization tensors are, by anal- 
ogy (see for instance fll~9l): 

ef ee ^(«-«) (13) 

The effect of such a polarized wave on a ring of test particles around a central particle 
is as follows: in the plane transverse to the GW-propagation,[] the ring is deformed into 
an ellipse that pulsates in and out along the x- and y-axes (e + ) or at an angle of 45° 
with respect to these axes (e x ). In Sec. |(| plane +-polarized GWs will be considered, 
entering a thin plasma. 

Unit circular-polarization tensors for gravitational waves are again constructed by 
analogy to electromagnetic polarization vectors, e# = (e x + ie y )/y/2 and = (e x . — 
ie y )/V2: 



e R = -j=(e + +ie x ) (14) 

e L = -y=(e + -ie x ) 

One can see that e^j and rotate a deformation of a ring of test particles in the 
counterclockwise and the clockwise direction respectively. 

An interesting characteristic of these polarizations is that the gravitational wave is 
invariant under a it rotation, which reflects the quantum behaviour of the associated 
massless graviton particles. 

The classical radiation field of a spin S particle is invariant under a 2ir/S rotation 
(for instance 2tt for spin-1 photons and 47T for spin-^ neutrinos). Apparently, gravitons 
are spin-2 particles. This also follows from the fact that the two orthogonal linear 
polarizations of a spin-S 1 field are inclined to each other at an angle \k/S. Indeed e + 
and e x are inclined to each other at jir (and for instance for photons e x and e y at ^7r). 



2.1.6 Stress-energy carried by GW 

A very important consequence of Einstein's equivalence principle is that local gravita- 
tional energy-momentum, or its density, is not uniquely defined. Non-vanishing stress- 
energy implies, after all, curved space, whereas the equivalence principle always allows 
one to choose a locally inertial (flat) coordinate system. Gravitational energy is therefore 
non-localizable. One has to average over several wavelengths to arrive at a well defined 
energy content in a certain region and talk about an effective smeared-out stress-energy 
of gravitational waves. 

In an arbitrary gauge, this effective stress-energy for a GW follows from contracting 
the weak field equations @ with h^ p leading to (see 35] p. 94 or [18| p. 955): 



1 in the reference frame of the central particle 
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„4 i 



which in the transverse-traceless gauge discussed before, reduces to: 



T ^ = -^G^ h Z) ( 16 ) 



This is the expression that will be used in estimating the amount of gravitational energy 
that can be converted into electromagnetic radiation or vice versa when a GW travels 
through a strong electromagnetic background field in for instance ( p7[) in Sec. ||| 

For a plane, linearly polarized GW of the form h^ u = 3? {(A + e + ^ u + A x e XfJiU )e lk ( z ~ tS) }, 
the only non- vanishing components are: 

= jm* = _ Tr; = _ Tt7 = A^flA, | 2 + \A X \ 2 ) (17) 

lb7rG 
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2.2 Tedrad formalism 

In describing physics on curved space-times, it is customary to work in the socalled 
tedrad formalism. This formalism is usually the most efficient in computing curvature 
and many other quantities. Also it provides the closest connection to measurement. 
After all, a measurement made by a (for instance accelerated) observer is always made 
in its locally inertial reference frame. For this reason, most calculations in this thesis 
will, at some point, involve the introduction of a certain tedrad system to describe what 
physics is going on. This section provides a short overview of the concepts that will be 
used in the following sections. 

2.2.1 Spinors and cosmological models 

In the General Theory of Relativity, the only equations of physical interest are covariant 
tensor equations. Half integer spin particles, however, do not transform like tensors, but 
like spinors, so the general 'recipe' of replacing all physical quantities by tensors, and 
using covariant derivatives everywhere instead of partial derivatives, does not work any 
more as soon as electrons or other spin ^ paricles are involved. It is therefore not possible 
to describe the effect of gravitation on the fields of electrons and most other charged 
particles (see 0]). Using non-coordinate methods, however, allows one to incorporate 
spinors into calculations on the same footing as tensors. 

Another area where the tedrad formalism is of great importance is in the study of 
cosmological models. A lot of work on this is done by G. F. R. Ellis (|l^1) and H. van 
Elst (in (l3j ) . They set up these models in a 1 + 3 covariant description where all tensor 
components are projected onto the average velocity vector and the directions orthogonal 
to it. The six propagation equations and six constraint equations resulting from this 
can only be completed by putting them into tedrad form, where one then has additional 
Ricci and Jacobi identities for the basis vectors .0 

Since the problems that will be discussed in this thesis involve a electron plasma in 
a cosmological background, it seems natural to adopt the tedrad formalism to describe 
these phenomena. 

2.2.2 Decomposition to locally inertial frame 

From the equivalence principle one is always free to choose a locally inertial coordinate 
system at every point P (local metric rj^), which is related to a general non-inertial 
system, by: 



2 The tedrad form also allows the set of equations to be put into a symmetric hyperbolic normal form 
which is easier to solve. 




(18) 



X*u,^ uVa/3 
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2.2 Tedrad formalism 



One can easily see any other choice of general coordinates x'^ would have led to the 
same expression, with: 

\' a = = ^L^l = ^ HQ) 

/i dx ^ dx'^dx" dx'» v K ' 

Apparently, A ^ does not transform like a single tensor, but as four covariant vector 
fields, hence the name tedrad or vierbein. Using this tedrad, any contravariant vector 
field or tensor can be decomposed along the locally inertial coordinate axes: 

_ \(p)p v (20) 

where the brackets denote which tedrad vector is contracted, and illustrate that by 
contracting the tensor fields with the tedrad vectors, they are decomposed into a set of 
scalars .0 This is precisely the reason why the treatment of spinors (of for instance Dirac 
electrons) is no longer different from that of tensors. Both can be put into an action as 
scalar s. 

In a more informal way of speaking, what is meant by the tedrad decomposition, is 
that at each point in space, physical quantities are projected onto artificial, orthogonal 
axis that 'stick out' of the curved space as it where. Hence the name non- coordinate 
frames. 



2.2.3 General covariance and Lorentz invariance 

To construct an action from the decomposed tensors, both general covariance and 
Lorentz invariance of the action must be satisfied. The scalar components are, after 
all, defined with respect to an arbitrarily chosen ONF, so the action should be invariant 
under a redefinition of this frame, viz under Lorentz transformations. 

F^\x) = A^(x)F^\x) (21) 
G Qw) = A^ a {x)A u p {x)G^\x) 

The same conditions apply to the tedrad vectors themselves. 

Now, to construct an appropriate action, one does not only need the fields, decom- 
posed as scalars, but also derivatives of the fields. Since al the fields are decomposed as 
scalars, the derivatives also have to be scalars. At the same time though they have to be 
invariant under Lorentz transformations. The derivative that satisfies these conditions 



is (see fTj, [|| or ||]) 



V f - \l)^- x{a p X {l)-^ X {S) ( 22 ) 



3 Later on these brackets will be omitted for brevity, assuming that it is clear to the reader what the 
different indices mean. 
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Summarizing: just as one obtains General Relativity by replacing the special rel- 
ativistic action and field equation by tensor equations, and the partial derivatives by 
covariant deriveves, one obtains a coordinate free description of General Relativity by 
decomposing all tensors or spinors into scalars (except for the tedrad vectors itself), and 
replacing all the partial derivatives with the covariant and Lorentz invariant derivative 




2.2.4 The connection 

In a particular set of coordinates, space-time is described by a metric given by g fiu (x) 
and a connection, that supplies the differential properties of space-time: the Christoffel 
symbols. In the non-coordinate formalism, space-time is is locally flat (metric rf LV ') with 
respect to a particular set of tedrad vectors, and its differential properties are given by 
a connection in the form of the socalled Ricci rotation coefficients. 
The specific form of these coefficients follows readily from (22). 



v A,A a = T 7 Q/3 A 7 & (23) 

r a _ \ 7 d \ P 

In tedrad components, covariant derivatives can thus be calculated completely analogous 
to the tensorial form: 

dT 

V a T /37 = A ( £ ) -^-r p^Tgy-T ^Tps (24) 
This derivative will be used in all the plasma calculations is Sec. ||. 
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3 Estimate of interaction efficiency 

According to General Relativity, gravitational waves and electromagnetic waves propa- 
gate with the same speed and, in vacuum, obey the same dispersion relation. If there is 
a linear dependence between the waves, they can therefore resonate and transfer energy 
([|l5|]). In this section propagation of light (F^ u ) in the presence of a strong magnetic 
background field (F^ v ) that is constant in space-time, is considered. 

The stress-energy tensor for the combined field (F'^ = F^ u + F^ u ) consists of: 

1. the square of the constant field term that does not generate gravitational waves, 

2. the square of the field of the electromagnetic wave which does not create gw's 
either, and 

3. the interference term which is proportional to the background field and produces 
an oscillating source term for gravitational waves. 

The field equations for the resonant term, which is the only relevant part of the stress- 
energy tensor, are given by: 



□<^„ = — j- 



8G 
T 



4 ' v 



a 3 



(25) 



a ± ua 

C 1 



where the second line results from choosing certain convenient coordinates, for instance 
such that: 



fVWFcp oc (E-E^-B-B^) 
E (0) = and B_LB (0) 



with 



(26) 



For an incoming lightwave along the x-axis (wavevector = k x = u/c), with no 
absorption or scattering (b(x) = b), try plane wave solutions with amplitudes normalized 
to unit energy density (|l6|) and /c" w = fc" IW = k a , e.g. gravitational waves produced along 
the x-axis. 



too 
Too 



16vrG 



< (V,o) 2 > : 



16vrG 



i-(E 2 + B 2 ) 



except for the overall 'conversion factors' a{x) and b: 



(27) 



F, 



bVAirf^e 



f v hv = 1 and Phi = 1 



(28) 
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For slowly varying amplitudes, neglecting terms quadratic in d/dx and of course k^k 1 * 1 
0, the field equations reduce to: 



k-r 



dx 2 



dx 



which is a differential equation that can be solved for a(x): 
da{x) 



dx 
a(x) 



ib 
ib 



c 1 
AG 




UC V f X F^ a (x')dx' + a(0) 



(29) 



(30) 
(31) 



One can already see from this that the energy of the generated GW grows with distance 
(a* a | x). In the figures below, the polarization of the generated GW is shown for two 
different polarizations for the incident emw. In Figure B the magnetic wave component 
is orthogonal to the magnetic background field, as in <M): B EM J_B(°)j_k EM , and the 
resulting GW is x -polarized. 







B(0) 














4 



Ikl 



y-axis 



Figure 2: Orthogonal magnetic fields. 
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3 ESTIMATE OF INTERACTION EFFICIENCY 



One could of course also choose the incoming emw to have a magnetic component 
parallel to the background field. A similar calculation shows that for B EM ||BW_Lk BM , the 
resulting GW is +-polarized. This is depicted in Figure 0. 



y-axis 



B (1) i 



B(0) 



Ikl 



Figure 3: Parallel magnetic fields. 



Linear combinations of these electromagnetic configurations lead to circularly polarized 
GWs (or any other polarization). 



Efficiency. To calculate the efficiency, assume: 

♦ a static, homogeneous magnetic background field, F^°\ 

♦ a typical length- or timescale for the interaction region, L = Tc, 

♦ unit convolutions of dimensionless amplitudes, 

♦ and no incoming gravitational waves, a(0) = 0. 

The energy transfer efficiency from emws to GWs is then given by: 



a 



a(x) 



1^/7(0)2 7-2 
C 4 



(32) 



This efficiency is the same for emw =4> GW and GW emw conversions since the relations 
are symmetric under time reversal. For a neutronstar binary or magnetar with a large 
surface magnetic field F^ ~ 10 16 Gauss and an interaction region (where one can speak 
of plane GWs) from R\ = 180 km to R2 = 500 km, and a dipolar decay of the magnetic 
field, something of the order of 10~ 8 of the energy could be converted, which still might 
be substantial considering the huge amounts of energy released in supernovas and binary 
mergers. 
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Fluctuations in the cbr. A very interesting example of GW to emw conversions, 
could have occured in the early universe. The observed fluctuations of in the cosmic 
background radiation might be explained by gravitational waves travelling throught the 
magnetic field that, according to some cosmological models existed throughout space 
at that time. If a fraction of the GW energy could be converted into emws, this could 
explain the ~ 10~ 5 relative anomalies in the gbr. This will be discussed in some more 
detail in Sec. |8[ 
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4 Exact calculation in vacuum. 

Motivated by the result of the preceding crude estimate, this section proceeds to the 
exact, analytical derivation of a conversion of photons into gravitons and back into 
photons again. 



4.1 Photons to gravitons 




Figure 4: Photon to graviton conversion. 



Consider the arrangement shown in Figure f| ( p6| and [0]), where an electromag- 
netic wave is partially converted into a gravitational wave. The coordinates are again 
chosen such that the incident electromagnetic wave propagates along the positive x- 
axes. Furthermore, the waves are taken to be linearly polarized plane waves: B z = 
E y = Aexpik(x — ct). For a completely general electromagnetic background field the 
combined Faraday tensor is: 



/ E x E y + Ae ik ( x ~ ct *> E z \ 

-E x B z + Ae^~^ -B y 

-E y - AeW*-**) -B z - Ae ik ( x ~ ct ) B x { ' 

\ ~E Z By —B x / 

As in Sec. [3|, the stress-energy tensor consists of two squared terms, that do not produce 
gravitational waves, and the resonant term of interest, which in this case is: 
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A_ 

4tt 



l —(Ey + B Z 

(Ey + B z ) 

—E x 
V ~B X 



~{Ey+B z ) 
(T>Z + Ey) 

B X 



E x B x ^ 

E x B x 

?y) ~{By + E z ) 

(By + E Z ) ~(B Z — Ey) J 



jk{x—ct) 



(34) 



satisfying the usual conservation condition T^- v = 0. So for the ten independent com- 
ponents of the stress-energy tensor, the Einstein equations (^) read: 



v ik(x—ct) 



(35) 



where the subscripts /, II and III refer to the regions as indicated in Figure ^[ The 
amplitude is defined by: 



□ 




= 


□ 




= C/e 


□c 




= 



4AG 



( ~(Ey+B Z 



\ 



'(Ey + B z 
E x 
B, 



~(Ey + B z ) 
(Ey + B z ) 

—E x 
—B T 



E x 

— Err. 



B x 

— Br, 



\ 



(B z - 

~(By 



Ey) 
-E z ) 



(By + E z ) 

-(B z — E y ) ) 



(36) 



The general plane wave solutions satisfying (35) are: 



I E iM x — c i) „ iklx—ct) i /~t n —ik(x+ci) 

C w I K7T. e + B^ u )e v ' + Cf„,,/,e v ; 



2ik 

C/^(H eife( "" ci) 



(37) 
(38) 
(39) 



where the arbitrary constants A^ u ^ , -B^) , , still have to be determined from 

the coordinate conditions (4>^ )U = as in (||)) and the continuity conditions at x = — L 
and x = 0: 



(Pl/i \x=-L 
<f>IIp\x=0 
*>IH U ,l\x=-L 
011^, iU=o 



?IIH \x=-L 
hlli,l\x=-L 

W^iU=o 



and 



(40) 
(41) 
(42) 
(43) 



Some tedious algebra results in the following constants: 
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4.1 Photons to gravitons 
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_ e 2ikL + 1 



e -2ikL _ 1 ^ 
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-2ikL 
-2ikL 



1 
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(44) 



D 



L 

2ik 



( 1 
1 

1 

V i 



1\ 



1 / 



In these equations the indices in brackets denote the constant factors of the correspond- 
ing index positions in C^. 

If the length of the device is set to a half-integer multiple of the wavelength, the 
amplitude of the retreating gravitational waves vanishes (for L = nX/2 — > exp(±2ifcL) = 
1). In this case gravitational waves are only generated in region III, which are linearly 
dependent on the length of the static background field. The complete expression for 
these waves in terms of the controlled parameters (background field, incident EM-waves 
and length of the device) becomes: 



L 



The last line follows from: = 4>u — ' 4>a with, from ([}(]), (p^ 1 = 0. As a result 



of this, ( |45| ) fully describes the metric of the GW entering the next part of the detection 
device (Sec. |4.2j ). 



4.1.1 Efficiency 

From the equation for the metric (["]). one can obtain an expression for the energy flux 
([l5|) carried by the generated gravitational wave: 

W = t m = ^ + (h\ 2 - hl 3 f) (46) 



(J^j (LA) 2 {(B y + E z f + (E y + B z f) 



Apparently, there is no coupling when a linearly polarized EMW with components B z = 
E y , propagates through a background field with E y = —B z and E z = —B y . Because of 
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symmetry arguments, the same result will hold for an incident linearly polarized emw 
with By = E z . For a pure magnetic field (i.e. F^ u — > B z ) this reduces to the estimated 
relation derived in the previous section, except for an overall factor of 1/167T, whereas 
for a pure magnetic or electric field aligned with the x-axis, the energy flux vanishes! 
Note the quadratic dependence of the energy flux on both the length of the static field 
and the amplitude of the incident EM-wave. 
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4.2 Gravitons to photons 



4.2 Gravitons to photons 




Figure 5: Graviton to photon conversion. 



Let us proceed now to the second part of the detection device, where the incident 
gravitational waves are reconverted to electromagnetic waves (Figure |). 4 The same 
symmetry relative to the x-axis as before is used to solve the Maxwell equations, 

Fot/3,y + F-ya,p + F/3y,a = and W^gF^)^ = = (47) 

in the oscillating metric perturbation h^ v . Because of the symmetry, only the derivatives 
with respect to x and ct contribute (to first order, these are just partial derivatives, since 
the connection is itself proportional to h^ v ): 



^23,1 


= -^23,0 


= 


b ,1 


— b ,0 


= 


^13,0 


+ -^30,1 


= 


P 31 
b ,1 


i p 30 
+ r ,0 


= 


-^21,0 


+ -^02,1 


= 


b T 


i p20 


= 



which results in the first field components -F23 = F 23 and F 01 = E% ' . In the specified 
metric, the contravariant components for F^ v are derived by: 



F"" = gT^Fap (49) 



4 After travelling through the first background field, the remaining light is blocked by a screen that 
does not influence the propagation of the Gws. 
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= {rf*^ + r^h" 1 * + h^rf p )F aP + 0{h 2 ) 

so for F 23 : 

F 23 = 4°) + (4°) - B^)h 30 + (4°) + E^)h 12 (50) 

Similarly, the contravariant components of the other covariant differential equations are 
given by (using the symmetries in the metric tensor): 

F 21,0 + F 02,l = {E iO)_ B ^ )h22l + iB (0) +E (0) )h23A (51) 

F i3,o + F so,i = (E (V-BP)h 23A + (BW+EW)h^ 
which together with leads to the wave equations: 



□F 21 = (52) 
OF 20 = (4 )- J Bf)/ l22 , 11 + (4°)+Ff)/ l2 3,n 



ae ik(x-ct) 



OF 31 



OF 30 = (E^-B^)h 23 , u + (B^+E^)h 33tU 
— p e ik(x-cb) 



with, 



a = ~k 2 {(4 0) - Bf))i 22 + (Bf + 4°>)fe,} (53) 
& = -k 2 {(E^-B^ 23 + (B^ + E^)^ 

In the same fashion, one can express ( |5l|) in ae lk ^ x ~ ct ^ / ik and j3e tk ^ x ~ ct ^ /ik. The general 
solutions of (|4^), (|5l|) and (|52|) are just in- and outgoing plane wave solutions plus a 
linear term expressing how the EMWs grow in amplitude while the GW loses energy. The 
electromagnetic field components are given by the solutions: 



B z = + ^e ik(x ~ Ct) + Ae ik ^- Ct ^ + Be~ ik{x+ct) (54) 

. Ey = _ E (0) + f™_™\ eik{x - ct) + Ae ik(x-ct) _ Be -ik(x+ct) 
V 2z/c 2/c / 

£> — ffW _i_ @^_ e ik(x-ct) _|_ £< e ik(x-ct) i jj e -ik(x+ct) 

y y 2ik 

-E z = -Ef + - A ^ e ^(- C t) + cre<*(*-«*) - De"*^) 
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4.2 Gravitons to photons 



where the arbitrary constants A, B, C and D still have to be determined from continuity 
conditions on the oscillating field (i = x,y,z): 

En\ x =-L = Ein\ x= -L and Bu\x=-l = Eiu\x=-l (55) 
Em\x=o = Eim\ x= Q and Bjji\ x=0 = Bim\ x=0 

4.2.1 Tedrad system 

First, however, one has to realize that a measurement of the field components in this 
curved metric only makes sense in a local Cartesian coordinate frame (onf) on the world 



line of an observer, as was discussed in Sec. 2.2. It is therefore customary to employ a 



non-coordinate frame for an observer at rest with respect to the GW, the tedrad A/ \, 
(§, ||], HH and H]). For these tedrad vectors to form a orthonormal frame, they 
have to obey the orthonormality condition: 

In physical terms, this is nothing more than a local coordinate transformation to an 
inertial system. 

In this case, the first components of and A/L can be chosen to be zero, and the 
last component of Aj^ to be zero to remove the arbitrarity of a rotation around the 
x-axis. The other components then follow from (56). 

A {0) = ( 1 + T^OO, 0,0,0) , A( )^ = (-/loo - lAo,^20,/i-3o) (57) 

*fi) = (^io,l - ^hu,-h 12 ,-hi 3 ) , A (1 )^ = (0,1 + -/in, 0,0) 

1 1 

Aj* 2) = (/l20,0, 1 - -h 2 2,0) , A( 2 ) /Lt = (0,/ii 2) l + -h 22 ,h 23 ) 

A {3) = (^30 , 0,-/123,1 - 2 ^33) > A(3 )/Lt = (0,/li 3> 0, I + -/133) 

The decomposition of the electromagnetic tensor in this onf is: 

E(a/3) = E^ ' X(a)fj.\l3)u (58) 

which for the electric field components results in: 

E x = E^-B^h 30 + B^h 20 (59) 
E y = (l-^(/ioo-M) 

_ ( aX _ a \ ik{x-ct) _ ik{x-ct) , g -ik(x+ct) 

\2ik 2k 2 ) 
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+ E^h 23 + B^h 30 - Bjphw + 4 0) /ii 2 
E z = E<® ( 1 - ^oo ■ 



f I3x 13 
\2ik ~ W 

+ E^h 13 + B^h 10 -Bi°h 20 



\ gik(x—ct) ^-j^ik(x-ct) _|_ j~)^—ik(x+ct) 



and for the magnetic field in: 



B x = BW+Ej®hn-EWh 13 (60) 
By = 4°) (1 + ^(^33 + ^11)) 



(3x ik^x—ct) _|_ fn ik(x-ct) , p. —ik(x+ct) 

+ 2ik + + 



# 2 = Bf) fi + lfa + hn)) 



ik(x-ct) _ a ik(x-ct) _ n -ik(x+ct) _ p>(0)u 

2ik v 
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4.2 Gravitons to photons 



4.2.2 Wave component solutions 

Now that the electromagnetic field components have been put in a non-coordinate basis, 
the previously mentioned continuity conditions are applied. The values for the arbitrary 
constants can be computed by requiring that in regions I and III only outgoing waves 
exist (Aj = Cj = Bui = Din = 0). This is done in Appendix |A[ The solutions of the 
electromagnetic waves in terms of these constants are given below for the three regions 
under consideration (see Figure |5|): 



Region I 

E y = -B z = - (V 2ifcL - l) e - ikix+ct) (61) 

x { + \ {Bf) ~ E y )){ ^° ~ *») + (^ 0) + 4 0) )^ + (4 0) 6o + 4 0) 6 2 )} 

E z = +B y = i (e~ 2ikL - 1) e ~ ik (x+ct) 

x {^2 - \(4 0) + ^ 0) )te>o - £ 33 ) - (4 0) 6o - 4 0) 6 3 )} 

Region II 

= 40) _ gfr + L) eik{x _ ct) jfr-kct) sinkx (62) 
y y 2ik 

x { 2^ + \ {Bf) ~ " i22) + ( ^ 0) + B y )] ^ + {B * ] ^ + ^ 12) } 

= (o) _ g(g + g) e tt(«-ct) _ gin kx 

2ik 

x {^2 + - 4 0) )(eoo - e 22 ) + (£i 0) + bW)&3 + (4 0) £ 3 o + 4 0) ei 2 )} 

2ifc 

x {2F - ^i 0) + 4 0) )(^o - e 33 ) - (4 0) 6o - 4 0) 6 3 )} 

^ = g (0) + P(x + L) e ik{x-ct) +e ifa-kct) sinkx 

X { 2F - + - fw) - - ^ )Cl3) } 



Region III 



E„ 



E-. 



= aL Q ik(x-ct) 

2ik 



(63) 



2ik 
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If the second part of the device is set up the same way as before, such that the length of 
this part is a half integer multiple of the wavelength (L = nir/k = ^n\), the retreating 
electromagnetic waves vanish, and what remains are just two linearly polarized electro- 
magnetic waves propagating in the positive x-direction. This is just what one expects, 
as it is the reverse process of the emw to GW conversion in the first part of the detector 
(in accordance to the remark above (|33|)). 



4.3 Efficiency 

It is possible to set up an arrangement where one generates an electromagnetic wave, 
which produces a gravitational wave. The remaining light is then screened off, while 
the gravitational wave is reconverted into an electromagnetic wave. A measurement of 
the thus produced electromagnetic intensity would proof the intermediate existence of 
a gravitational wave. This intensity, in the magnetostatic case, is given by: 



W 



( GLV 
4^ V c 4 

4tt 



) [{B z Bf)+B y B^f + {B y B^-B z B^?\ (64) 
{LB z f 



where in the last line the length and the electromagnetic field of the two regions are 
taken to be equal. One can see that the final intensity depends on the fourth power of 
both the length and the magnetic field used (and on the square of the intensity of the 
initial electromagnetic field A). Still, the smallness of the factor (G/c 4 ) 2 does not make 
this a very feasible terrestial detection device at present. 

However, it will be shown in the next few sections, that the astrophysical relevance 
of the discussed GW to emw conversion is significant. In most cases, GW to emw 
conversions will occur (and not the reverse process) the efficiency of which is given by 
the square root of (|64|), For a magnetostatic background field, this, of course, results in 
the same efficiency as in Sec. |3| and Sec. |4.1. 1[ 

The main purpose of the calculations in this section was to rigorously derive the 
characteristics of GW 44> emw conversions in the most general circumstances. 
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5 Magnetars and binary mergers 

In the previous sections, the possibility of converting gravitational wave energy into 
light in the presence of a static, homogeneous electromagnetic field was discussed. Such 
a conversion is interesting for two reasons. First, it could be relevant for the energy 
excess in the electromagnetic spectrum of very energetic phenomena such as gamma ray 
bursts. The reason for this is that a sub-class of GRBs is supposed to be powered by 
binary mergers, which release most of their energy in GW and not in EM radiation. 

Secondly, it could be a useful indirect detection device of gravitational waves, because 
emws are much easier to detect then GWs. 

One should therefore look for astrophysical sources with the following ingredients: 

1. Obviously, one would like to have a source that produces gravitational waves that 
carry large amounts of energy, so that even if only a small fraction is transferred 
to EMWs, the result is significant, 

2. The source has to produce GWs with relatively high frequencies (~ 10 kHz). Oth- 
erwise, the emws have a frequency (the same as the GW) below the interstellar 
plasma frequency and will be absorbed before reaching earth, 

3. The interaction will only take place in either an extremely strong magnetic field, 
or a (weaker) field that prevails over extremely large distances (i.e. a primordial 
field). 

4. The region of interaction is a vacuum or a thin plasma, so that one can neglect the 
difference of the dispersion relation for the emw from vacuum and GW-dissipations. 

The most promising candidate, satisfying these demands, is the so called magnetar, a 
compact neutron star formed by an a-symmetic supernova collapse, a merger or collision 
of two neutron stars as studied by Kokkotas et. al. (pl|-pO[|). 

What remains after such a collapse, is a violently oscillating neutron star, quaking 
in several distinct quasi-normal modes. Post-Newtonian fluid modes, such as the f- 
mode (fundamental mode), p-mode (pressure as restoring force) and g- modes (gravity 
as restoring force), are dominated by the fundamental mode with typical frequencies of 
1 — 2Hz. The socalled u>ave modes, on the other hand, are purely general relativistic 
modes. These are not fluid modes, but GWs due to strong space-time dynamics, with 
typical frequencies of 8 — 12kHz (and up to 38kHz). This means that the w modes 
might be able to produce EMWs in the right frequency domain (even without photon 
acceleration) . 

The energy radiated by these GWs is a fraction of a solar rest mass energy and 
is radiated in a small damping time in axisymmetric waves, producing a very large 
gravitational flux density. 

Finally, the typical magnetic fields of such magnetars are the largest ever encoun- 
tered: ~ 10 16 Gauss. If one now assumes that the far field of the magnetar can be 
regarded as a vacuum (Sec. 5.1) or a thin plasma (Sec. |6|), all our conditions are 
satisfied. 
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The next subsection gives an estimate of the produced electromagnetic energy flux 
due to gravitational waves in the presence of a magnetar using the results of Sec. [|| 
and Sec. 4.2. Then in Sec. || the same will be done for a thin plasma using a slightly 
different approach in which the Maxwell equations are already put in a non-coordinate 
basis before solving them. This allows easy comparison of the current densities in the 
plasma with the gravitational effects induced by the GW. 
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5.1 Estimate of efficiency for magnetar 

In Figure ||, the (far) magnetic field of a magnetar with mass M = Mq, radius R\ 
Rs = 2GM/c 2 = 3 km and surface field B\ = 10 16 Gauss is shown. 
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Figure 6: Gravitons to photons around magnetar. 



The magnetic field of such a neutron star will be that of a dipole out to the light 
cylinder and a spherical 1/r decay beyond that. Since, however, the studies [31]-[f2C|] 
involve a non-rotating star, only the dipole field will be of importance. 

Furthermore, one can only speak of a GW in the radiation zone, for r ^> R±. In this 



estimate, following Brodin, Markhmd et. al. [11], it will be assumed that the interaction 



becomes effective around i?2 = GORs = 180km. 

B(r) = B 1 (—) , Rx<r<R 3 (65) 



To compute the amplitudes of the generated EMWs, use (|32| ) or fl46|) . To apply this 
formula to the reverse case, realize that the efficiency of GW to emw conversion in a 
static magnetic field is the same as for the reverse case, since the relation connecting 
<fifiv ~ and F^ v is linear and symmetric under time reversal. So, for c = G = 1, 
b = E = —B, and a(z)/k = h, where h is the dimensionless effective amplitude of a 
+-polarised plane GW, the efficiency is: 



h 1 2 

And using the GW amplitude as calculated by Kokkotas et. al. in JT|: 
h m-2i ( E V /2 ^lOkHzV/ 2 /50kpc\ 

the EMW amplitudes can be calculated: 

E y (r,t) = -B x {r,t) (68) 
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The most energetic magnetars, with E = 10 -2 M@c 2 will have incident GWs at r = 180km 
with an effective amplitude of h = 0.001. Putting this into the equation for the produced 
electromagnetic waves and choosing some cut-off distance i?3 where the magnetic field, 
or just the coupling, becomes negligible (or 1/-R 2 is just negligible with respect to 1/-R 2 ), 
the result is: 

E ymax = -B xmax ~ 10 12 V/m (69) 

using k = to with lo/2tt = 10kHz. 

This large value results mainly because of the extremely large magnetic field, which 
is four orders of magnitude larger than for instance the field of a neutron star binary. 



5.2 Comparison with neutron star binary 



If, in (66), a close binary just before merger is considered as a source, with the same 
distance scales, but a magnetic surface field of 10 12 Gauss, an effective (Newtonian) GW 
amplitude of 4 • 10 -4 and a GW frequency of 2 • 10 3 rad/s (caused by the rotation of the 
binary instead of the magnetar star-quakes), the emw amplitude is still very significant 
with 50MV/m. 

The frequency of the latter emw, however, causes these waves to be damped by 
the interstellar plasma, whereas the former could easily be detected by a space based 
detector (avoiding the ionospheric cutoff), as long as the source resides close to our 
galaxy. In fact, for the extreme case discussed here (maximum magnetic field and 
maximum gravitational wave energy), the emw amplitude reaching earth from a source 
at a distance of ~ 50kpc would be as much as lV/m if the emw only decayed by spherical 
attenuation. 

It has to be realized, though, that the values mentioned here are probably highly 
exaggerated when one takes into account all the real dynamics around a magnetar or 
neutronstar binary. More research has to be done to deal with spherically decaying GWs, 
non-uniform plasmas etc. 



5.3 Improvements to calculation 

As a first step towards a more realistic astrophysical setting, the vacuum has to be 
replaced by a thin plasma. In the next sections ( |6.1| -|6.4|), it becomes apparent, that 
the interaction is hardly effected by the presence of such a plasma. The presence of a 
plasma induces a slight wavenumber offset as compared to the vacuum waves, but as 
long as the coherence lengths of the GW and emw are longer than the extension of the 
static magnetic field, the equations for E y (r,t) = —B x (r,t) will remain the same as for 
the vacuum case. The main importance of adding a plasma is to either damp or excite 
the EMWs to higher frequencies. 

In Sec. 0, the plasma will be considered in the magnetohydro dynamic approximation 
of a perfectly conducting ideal fluid with pressure gradients. A dispersion relation is 



derived in Sec. 7.4, in search of plasma waves such as Alfven waves, magnetosonic waves 



and slow and fast modes, excited by the GWs. 
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6 Plasma calculations 

To obtain the expressions for the EMW produced by a GW traveling through a plasma 
in more or less the same fashion as in Section [|, the Maxwell equations are written in 



terms of a general local orthonormal frame. This will be done in the next section S.l 



In Sec. S.2 the ONF will be specified for a plane polarized GW. 



6.1 Non-coordinate Maxwell equations 

The general, coordinate, decomposition of the Faraday tensor with respect to an observer 
moving with a velocity u a , satisfying u a u a = —1 is: 5 



F ab = u a E b _ u b E a + e abc Bc QT ( ? q) 

E a = Fo a , B a = -e abc F bc 

where e abc = u d r] abcd is the volume element on hypersurfaces orthogonal to u a , which is 
skew symmetric in all its indices and satisfies £123 = f?oi23 = y/\g\ = 1- 

To evaluate the Maxwell equations in a locally inertial non-coordinate system, again 
an orthonormal tedrad is used: eo = u and e a = V a (with e^e b ° gij = r] a b)- I n this ONF 
the covariant Faraday tensor is just: 

F ab = u a E b - u b E a + e abc B c (71) 

and the covariant derivatives are defined, completely analogous to the usual tensor form 
in terms of Christoffel symbols, as: 



V a F bc = e a (F bc )-T d ba F dc -T d ca F bd (72) 

da^ 



V a F bc = e a (F bc )+T b da F dc + T 



where in the non-coordinate desciption, the Christoffel symbols are replaced by the 
socalled Ricci rotation coefficients. These are defined by: 

r c a6 = eW V J e a Ve 6 e a = r c afe e c and T {ab)c = (73) 

The general form of the inhomogeneous Maxwell equations in the ONF can now be 
expressed as: 



V a F ba = e a (F ba ) + T b da F da + T a da F bd (74) 
= e a (u b E a )-e a (u a E b ) + e bac e a (B c ) 



5 From now c = 1, roman indices are 0,1,2,3 and greek indices 1,2,3 following the conventions in 
0. 
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+ T b da u d E a - T b da u a E d + e dac T b da B c 
+ T a da u b E d -T a da u d E b + e bd °T a da B g 

and similarly for the homogeneous equations. In evaluating these expressions, one has 
to realize that: 

1. E a = E a and B a = B a , so all terms in E° or B° vanish, 

2. the time axis is chosen such that eo = u, so the observer is at rest with respect to 
the ONF and only u° = 1 is nonvanishing (all terms u a are zero), 

3. e 0ab = u d rf abd = u V 0ab0 = and similarly e a0b = e ab0 = 0, 

4. since the rotation coefficients are skew symmetric in the first two indices T aa f, = 0. 

Inhomogeneous equations Using these 'selection rules', one finds for the first inho- 
mogeneous Maxwell equation (V a -F 0a = 4-7rp): 



V a F 0a = e a (E a )-T so E s + r so E 5 + T a Sa E 5 + e 5 ^T Sa B 7 (75) 

v-e = (-rv^-e Q/37 r% J B 7 ) + 4vr Pm 

= pe + 47rp m 

where the last equation defines the extra 'charge density' pe- The other three inhomo- 
geneous equations follow from V a F^ a = 4k 



V a F? a = -eo(^) + e^e a (B 7 ) (76) 

+ (rV^-r^-rv^ 

+ (e Sa ^ Sa + e^T a Sa )B, + e^T° S0 B, 



eo (E)-vxB = (r a 0/3 -r%)^-r^ 

+ e^(T% B 1 + T 5 ^Bs)-^i m 
= -Je - 47rj m 

where the intermediate result that {^ Scr/ T g a + ^T^jB^ = e^TVS a is proved in 
Appendix M, and where the last equation defines je- 
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6.1 Non-coordinate Maxwell equations 



Homogeneous equations The first homogeneous Maxwell equation follows from 
chosing (abc) = (a/37): 



V(aF /37 ) = Vai^r + V 7 F a/J + VjgF 7a (77) 
= ef3yS^a(B S ) + e a /3se-y(B s ) + e ya se f3 (B s ) 



■pd 77 77 "pd 771 



-es(B s ) - T° a/3 E y + r° 7Q .Eg + Y ^E a 

- 1 - 1 — 1 7/ 0- c 'a 

+ e 7<5r>r <5 /3Q ,fi I/ + eps v T s ai B u + e a s u T s 1/3 B u 

— £aSuY S /3ryB u — e^SuY S a fjB u — eps v T & ^ a B u 
= e^r%F 7 - e^e^r^ 

= e ^r%F 7 + r%ir (+r^ = o) 

v b = -rVB^^r> 7 = Pi? 

The other homogeneous Maxwell equations follow from chosing (abc) = (0aj3): 



-e (B) + V x E 



e (B) - V x E 



(ep(E a ) - e a (Ep)) + e aPl e (B^) 
(X d a0 F d/ 3 + T d p F a( i) 

(^ d oi3 F da + r a/3 F od ) - (r^Fdo + r d 0a F^) 

"-( r °a0-^/3 ~ r %Fa) - r nfl F n 



(r 5 a/3 F 5 



/9a 



F 5 ) 



0/3 -^a 
r °0aF/3 



+ ( e / 35 7 r <5 a0 - eQ^r^o)^ 7 



e^(T% E y + r 5 pi E 5 ) 

e I//3a e Q , ( 5 7 r 5 o / 3)-B' y 
<5 



t l e e^T 1 0/3 



o/3 7 / r 
-j-B 



f 7 + r 01 e§) 



p/3 R a 
1 0/3"° 



(78) 
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6.2 Gravity induced charge and current densities 

In this section, the ONF used in the non-coordinate Maxwell equations (f75"l)-(f78|), is 
specified to a plane, +- polarized GW with metric: diag(— 1, (1 + h), (1 — h), 1) (compare 
with (12)). The natural tedrad for this system is, as one can easily verify by applying 
the orthonormality condition ((56|): 



e (o) 



(1,0,0,0)) 

(0,(1-^/0,0,0) 

(0,0,(1 + ^,0) 
(0,0,0,1) 



(79) 



In this ONF the only nonzero rotation coeffiecients are T 1 



13 



-F 2 — LQh. nnrl T 1 
1 23 — 2 dz allCL 1 10 



20 



^^r, since only the z and t derivatives of the e2 and e3 vectors contribute. 



This imediately tells one that the gravity induced charge densities vanish. The only 
components that remain for the induced current densities are: 

j E = [T^ E^-e a ^T s ^Bs]e a (80) 

j B = [r a po BP - cooe a ^T s M E 5 ] e a 

Explicitely, the second term in the equation for for a = 1, (3 = 2 and 7 = 3 becomes^: 



(r s 23 -r 5 32 )B s = r s 23 B s = (o,( 



( (1 + U)5 a 
B z 



1 



hB« 



(81) 



whereas for a = 2, /3 = 3 and 7 = 1: 
( p5 3i ~ ^ S 13 )B 5 = 



T 5 13 B S = -((1-1^,0,0)^ 



;i + \h)B x 



- - 2 h)B y 
B 7 



= --hB x 
2 

Similarly from the equation for j# (replacing B — > — E) it follows that: 



T S 23 E S 
T 5 3i^ 



—hE v and 
2 y 

\hE x 



using h 



dh 

' at ■ 



(82) 



(83) 



JOACHIM MOORTGAT 



WATCHING GRAVITATIONAL WAVES 



37 



6 PLASMA CALCULATIONS 



6.2 Gravity induced charge and current densities 



and for the other components of j#: 



and for j^: 



2dt x 
ldh F 

2dt y 



(84) 



r 1 ^ 



2dt 
l9h B 

2dt y 



(85) 



Putting things together: 



1 



3e 



3 b 



By - E X 
B X + Ey 






(86) 



The Maxwell equations in the presence of a gravitational wave are summarized below 
in their final form for future reference. 





V • E 


= 47rp m 


(poisson) 


dE 

at 


- V x B 


= -47rj m - j E 


(ampere) 




VB 


= 


(no monopoles) 


dB 

at 


+ V x E 


= ~3b 


(faraday) 



With }e and defined in (Bfj). 
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6.3 Energy-momentum conservation 

To evaluate the Maxwell equations, or to find electromagnetic wave equations, one still 
needs an expression for the matter current. This expression follows from the linearized 
equation of motion. In this case, for a one component plasma where the ion veloc- 
ity is negligible compared to the electron velocity, one can use the energy momentum 



conservation equation (see |L6| or p3[). 

The matter part of the stress energy tensor, for a pressure free 'dust-like' plasma, is 
(see Q and |2§) Tf = pV a V h (where V a = (1, v) is the non-relativistic fluid velocity 
for 7 <C 1). The electromagnetic part is not needed here because of the useful identity 
V^T^f = —F ab (j M )i). Thus the conservation equations follow from: 

V b {mnV a V b ) - qn(u a E b - u b E a + e abc B c )V b = (87) 

or, in a 1 + 3 split, 



|^+V-(nv) = -n(T a 0a + T% v Q + T^y: 



^ + v-v) v = i(E + vxB)-(r 00 + (rv + r%K + r>V) 

In the specified metric, these equations can be linearized around the unperturbed plasma 
state: E = 0, v° = and v(z,i) = v(z)e _tw *, where all the barred quantities are first 
order amplitudes. In that case, all the terms containing the rotation coefficients are 



either zero or quadratic in h, which reduces the equation of motion (88) to: 

<9v o - 

— = —icov = — (E + v x B) (89) 
at m 

Now take j m = qnQV (with no the unperturbed electron number density) and intro- 
duce the following relevant frequencies: 

CYCLOTRON FREQUENCY: =^ UJ C = (90) 



PLASMA FREQUENCY: 
UPPER HYBRID FREQUENCY: 
and solve (^9|) for v or j M : 



LO c = 


qB 

m 


UJ p = 


/ Aimoq 2 


V m 


= 





1 ui, . 

(jm)x = -. ~iE x (91) 

47T LO 

&0, = r-K + - J *)(^V) 

47T LO V CO / \LO z — L0£ / 

47T LO V LO J \LO z — LO Z c ) 
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6.4 Exact calculation for the plasma waves 



The component needed in the next section is: 



4muj(j u ) y = -ui 



U3 n 



En 



(92) 



6.4 Exact calculation for the plasma waves 

Now that the Maxwell equations are put in an elegant and transparent form, it is 
straightforward to obtain wave equations for an EMW excited by a GW passing through 
a static magnetic background field B^ ) = (B(z), 0, 0). For the electric wave components, 
take the curl of Faraday's law ([Tq ) and substitute it into Ampere's law, differentiated 
with respect to time. And to find magnetic components, proceed the other way around. 
Thus, one finds: 



(J^-V*\e = -^(4jrj„+j E )-VxjB (93) 

£- v2 ) B = «"*.+« -t!f 

To solve these equations assume the GW to be monochromatic and let the frequency of 
all perturbations coincide with that of the driver, leaving the spatial dependence free 
for the moment, i.e.: 

h = he^-V (94) 
E = E(z)e~ iujt 
3m — 3m{z)e 

B = B(°)(z)e 1 + (B)( 1 )(z)e-^ 



Dropping all second order terms in the perturbation (which are in effect just the current 
density terms in h • E), the relevant wave components are: 

(J^ + co 2 ^E y + 47riu(j m ) y = -^k 2 B-±ik^Jhe ikz = F(z) (95) 
d 2 \ 8 ( 1 8b\ 

9^+^)B x +^-(jm)y = [k 2 B--ik—)-^ = -F(z) (96) 



where F(z) is a function of z through the specific form of the static background field. 

The expression for E y is now easy to evaluate using (p2), while B x is harder to find, 
since the z-dependence of (j m ) y has not been specified. B x can, however, be deduced 
from E y , as will be done below. 
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Solution for E y : Inserting ( |92D in (p5|), leads to an inhomogeneous differential equa- 
tion for E: 

/ d 2 



dz 2 



i 2 2 
+ UT - LOp 



U 2 ~U 2 

oj 2 - CJ? 



E y = F(z) 



(97) 



To find the solution of this equation, one first determines the solution of the homogeneous 
equation, taking F(z) = 0, and its Wronskian: 



E H y{z) 

W 



Ae tklZ + Be 
Ae iklZ 
k 



(98) 



det 



ik,Ae iklZ 



Be~ iklZ 
-ikiBe~ iklZ 



2ikAB 



In (|98|), ki is the wavenumber of the emw satisfying the dispersion equation kf 



(Aw) which is slightly shifted compared to the dispersion relation in vacuum. 



This shift is given by: 



Alo 



UJ 2 — U! 2 



(99) 



The general solution of the inhomogeneous equation is then given by: 



Ey( Z ) 



2k 



iz 



-ikiz' 



lZ F(z')dz' + e~ lklZ / e iklZ F(z')dz' 



o 



(100) 



Here, as in Sec. ||, the assumptions have been made that the static background field is 
restricted to a region < z < L and that there are no incoming EMW's in the interaction 
region. 

From ( |100P , the amplitude of the outgoing wave can be calculated by inserting the 
expression for F(z), expanding terms in Alo <C oj and assuming that k\ is approximately 
constant in the interaction region. 



(E 



yjout 



2 - 1. -iAoj 2 

< oj h %0Jh 

2k, 1 2 2oj 



B(z)e\ ^ 



'dz 



(101) 



ihoj \fk~i 



2 h Jo 

ihoj f L «, N 



'dz 



'dz with 



h 



(102) 



where the second term in ( |101|) follows from partial integration of and the stockterm 
vanishes at the boundaries because the magnetic field is localized. The full solution for 
E y (z, t) is then just: 



E y (z,t) = (E y ) mt e** 



i(klZ — LUt) 



(103) 
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6.4 Exact calculation for the plasma waves 



Solution for B x : The calculation of the magnetic part of the outgoing wave follows 
realizing that all that: 

4n-j^(jm)y = 4niki(j m ) y ^ Amu ( 1 - i (~^) ) Um)y ~ ^■niu}(j m )y (104) 



Comparing this with ( |95D and |96| ), it follows that B x ~ — E y , as one expects. 

The emws, or 'light', generated by the gravitational waves are apparently hardly 
effected by the presence of a thin plasma. The wavenumber has a small offset with 
respect to the vacuum solutions, reflecting some dispersion caused by the plasma, but due 
to the strong magnetic background field, this will be a completely negligible effect. The 
wavenumber 'mismatch' is after all inversely proportional to the cyclotron frequency (in 
the regime u>p/uj c <C uj ~ lo p <C u> c ) and thus also to the magnetic field, so for a magnetic 
field of say 10 8 Gauss and a typical frequency of a few kHz (Aw) 2 /uj ~ 10~ 19 rad/s. 
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6.5 Magnetars revisited 

As mentioned in the previous section, the presence of a plasma in the interaction region is 
not very important. This is the result of the fact that the large magnetic field suppresses 
the mobility (resulting in Au oc 1/uj c °t *C 1). Assuming therefore that the 

coherence length of the sc gw and emw is larger than the length of the interaction 
region, the same results are found as for the vacuum case. The difference lies in the 
fact that in order to travel over astrophysical distances, the generated emws have to 
overcome the interstellar plasma. As soon as the magnetic field has decreased such that 
uj c < u>, (|89|) and a spherical decay of the EMW beyond the interaction region allow the 
electron quiver velocity to become relativistic at a distance: 



v(r) ?-—E max (105) 

rmo r 

r rol ~ 3 • 10 12 km 

and the electron velocity somewhere in between. In that case the EMWs become highly 
nonlinear and effects such as parametric excitation and, even more importantly, har- 
monic generation become important, which is interesting, because the latter effect can 
convert the emws to even larger frequencies. 

It is not likely that this, non-linear, effect will play a very important role in the 
interstellar plasma, exciting higher harmonics with 60 times the frequency of the original 
waves, as has been done in laser experiments in a laboratory plasma. If, however, only 
the first or second higher harmonic is reached (with the same energy), this would be 
enough to overcome the interstellar plasma damping, and the resulting light should 
be detectable with a telescope such as the proposed Astronomical Low Frequency Array 
(alfa). Supernova? and collapsing binaries within the Local Group and the Virgo cluster 
might be detectable in this way with a event rate of as many as a few per year as soon 
as detectors such alfa are operational. Furthermore, these signals could be compared 
with the GW observations when the first GW detectors such as ligo, LISA and VIRGO 
open the GW observation chanel. 
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7 Magnetohydrodynamics 

Another way to study interactions in a plasma is of course in the magnetohydrodynamic 
approximation. Consider an infinitely conductive hydrogen-like plasma in a singe fluid 
approximation (for more details, see jn| and |23[]) where a pressure term is added in 
comparison to Sec. |6.3| . In other words, consider the full energy-momentum tensor. 

In an incompressible plasma, the unperturbed pressure will be negligible with respect 
to the energy density fi = p(l + II) (with pH the, small, internal energy density), but the 
pressure gradients will be important. To avoid mistakes, all terms will be maintained 
throughout the exact calculation. 



7.1 Conservation equations 



The stress-energy tensor for a plasma in an electromagnetic field can naturally be de- 
composed in the tensor for the energy in the pure electromagnetic field that of the matter 



field, as in Sec. 3. J 



-tab 



rpab _j_ rpab 
J m ~r em 

(ji + p)V a V b + pr] ab + -L 

47T 



(106) 



F ac F c b - -r] ao F ca F cd 



ab jpcd ] 



which makes the conservation equations easy to evaluate (using the result of Exercise 
3.18 in M): 



v T ab 



V 6 T" 6 + V,T" fe 



V fe T r 



ab 



F ab j b 







(107) 
(108) 



Resulting, using the 3 + 1 split, in the equation of energy conservation (a = 0) and 
the equations of motion (a = a = 1, 2, 3). Again V a = (1, v) T and j a = (p m , jm) T - 



d dv 



jE 

(p + p) 
P 



r° db v d v b + v b db v d 







o 



v ( Ql^ + P) +V ' (^+P) V ) + (P+P) ( ^7+ v - V ) v + Vp 



Ot 



(109) 



(110) 



j m x B - pE - (v • j m )E + (j • E)v 

(p + p) 



r a db v d v b + r b db v a v d 



V 



r dbV + r db v 
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Most terms in these rather long expressions vanish in a first order calculation in the 
metric of a plane polarised GW (ffs|): 

o All the rotation coefficients in the first equation vanish by themselves in this metric, 

o the only non-vanishing Ricci coefficients in the second equation result in: 



-i(fi + p)(k g v z + uj g )h | -v y ) (111) 



2 



which is clearly non-linear in the perturbed quantities, 

o both E and v and j m have no zeroth order components in the perturbation, so all 
the products of these quantities are second order or higher, 

o pEi is negligible with respect to j m x B due to the charge neutrality in a plasma 
(the global charge density is very small, whereas the currents can be substantial; 
for details, again, see [16] or p3|), 



o for the same reason the convective term v • Vv vanishes. 

A final simplification is made by substituting the first equation in the second. What 
remains are the following simple conservation equations: 

^ + (/i + p)V-v = (112) 
(/* + P)^T + Vp = j m xB (113) 



dt 



-L(V x B) x B - ±-] E x B (114) 



where in the last line, 'Amperes' law is inserted without the displacement current, which 
is negligible with respect to the current j m ~VxB^>^-^in the MHD approximation. 
The above equations are just the usual mass conservation equation and the equations 
of motion, but now with added terms that express the gravitational effects. 

7.2 Ohms law 

To close the system of equations, one also needs the generalized law of Ohm, which 
follows from the separate equation of motion for the electrons in the approximation that 
the electron inertia is negligible {dv e /dt = 0). This is adequate as long as the pertur- 
bations caused by the GW are slow compared to the motion along the magnetic field (a 
weak constraint for the cyclotron motion around an extremely strong magnetic field) , so 
the electrons have ample time to reach dynamical equilibrium. For zero resistivity, the 
single- fluid electron equation of motion is then given by (compare to ( |110| ) and ( |114[) ): 

Vp e = /) e (E + v e xB) (115) 
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or, from j m « p e (vj - v e ) and therefore v e « v - v* ps v - j/p e : 

E = —v xB + — (jxB — Vp e ) (116) 

However, the more trivial and commonly known, Ohms law E = — v x B prevails in 
the limit of a small Larmor radius with respect to the typical lenght scale of the GW 
perturbations (|p3fl). Since the magnetic field is very strong, this Larmor radius will 
be smaller than 10 _8 cm for electrons moving with thermal velocity {v^ e ~ 6 • 10 7 m/s 
at T ~ 10 8 K) and still smaller for lower perpendicular velocities. Even the smallest 
expectable GW perturbations should therefore be much larger than the Larmor radius. 

7.3 mhd Maxwell equations 

The reason for all this is to eliminate the electric field entirely from our set of equations, 
which is achieved by inserting the linearized, general law of Ohm in Faraday's equation: 

<9B 

— = Vx (vxB)-Jb (117) 

Linear gravity terms The gravity induced current densities used in this equation 
and in ( |114| ) are to first order: 




(118) 
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7.4 mhd wave dispersion relation 

Finally, it is possible to derive a dispersion relation for the plasma waves, generated by 
the GWs that perturb the plasma. To obtain such an equation for the combined matter 
plus electrodynamic energy density content, the Faraday tensor no longer suffices, as 
it did in Sec. [|. What is needed now is the second covariant derivative of the total 
stress-energy tensor, which amounts to substituting (minus) the time derivative of the 
energy conservation equation (|113| ) into the divergence of the equations of motion (|114| ) . 
This is somewhat similar to the procedure used when deriving magnetic wave equations 
by substituting the time derivative of Faraday's law into the curl of Amperes law (or 
the other way around to obtain electric waves) . 
To first order, the wave equation becomes: 



V a V b T ab = — ^ (v b T ob ) + V- (v b T ab ) (119) 



dt 

^-V 2 P = -^•(VxB)xB + iv.(j £ xB) 
= ^£ V 2 B-^V 2 h 

47T 87T 



or, with the definition of the sound velocity: 



To eliminate B, differentiate ( |117|) with respect to time and insert V-v from the equation 
of mass conservation. 



d 2 B ^d^ __ dj B 



B^V-v-^f (121) 



dt 2 dt _ dt 

= b£/5 + 5 



dt 2 \x 2 

Before continuing, first define the relativistic Alfven speed by: 

1 1 / B 2 V 1 Bl u\ , . , 

^ = ? + (4^j °' 4^ = T^t7? W, ' h X = " + " <122> 

This expression reduces to = u A as long as we are discussing the mhd approximation, 
where the displacement current is negligible with respect to the current density. One 
has to keep in mind though, that in the limit of a very thin plasma (j vacuum) with a 
strong magnetic field, the asymptotic behaviour of u 2 A is that it goes to c and of course 
not to infinite speed. 

To arrive at the final dispersion relation for plasma waves excited by a incident 
gravitational wave, allow for damping, so the energy of the GW can dissipate into the 
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plasma. In other words, keep the frequency and wavenumber of the plasma waves 
unrestricted and try perturbations of the form: 



h = hS kaZ ^ st) 

B tot = B + Be^^-^ 



(123) 



Inserting this into (|120|) and (p_21[), results in the dispersion relations: 



{-u 2 + k 2 (c 2 + u A )} fl = ^ 



k 2 



LOg U! 2 



h 



(124) 



A more physical form of this equation is found after the left-hand-side of the equation is 
made dimensionless (just as the right-hand-side in h) by normalizing the energy density 
perturbation to the total energy density: e = jl/x. 

Dividing both sides of (124) by x, leads to the final dispersion relation: 



{-^ + k\4 + u A )}e 



Ui 2 in 2 

2 U A LO g 



r *2 k*\ 



h 



The characteristics of this equation are discussed in the next subsection. 



7.5 Summary 

In this section, the equations of total energy density conservation and motion were 
derived in a tedrad system describing the influence of a passing gravitational wave on 
the underlying metric. It was found that these equations, to first order in h, have exactly 
the same form as in a flat metric. 

The coupling with the GW comes in through Maxwell's equations, used to eliminate 
B, with extra gravity induced terms. Standard MHD approximations of infinite conduc- 
tivity and negligible displacement currents were used to eliminate E from the equations. 
Finally, a wave equation was derived from the second covariant derivative of the total 
energy-momentum tensor for a perfect fluid in a electromagnetic field. 

The dispersion equation that resulted from this, satisfies the expected limiting be- 
haviour: 

♦ In the limit of vanishing GWs, h — > 0, it reduces to a common fast magnetoacoustic 
plasma wave with lo 2 = k 2 (c 2 + u 2 A ). 

♦ Without damping (when k = k g and u) = uj g ), the interaction disappears, and one 
finds the same magnetosonic solution. 

♦ In the limit of vanishing density, e — > 0, the displacement current is no longer 
neglegible, and one has to use the relativistic expressions for 11,4 and c s , that go to 
c and respectively in the vacuum limit. The 'plasma' dispersion relation therefore 
tends to k = uo (vacuum solution) and one retrieves the dispersion relation for a 
gravitational wave in vacuum, k g = uj g . 
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♦ Finally, the strength of the coupling depends on the square of the magnetic back- 
ground field through the Alfven speed and on the square of the frequency of the 
driving GWs (compare to the linear dependence on the frequency of the amplitudes 
of the emws generated by GWs). 

For a class of gamma ray bursts, powered by merging neutron star binaries, it would 
be very interesting if even a small fraction of the enormous amounts of energy released 
in gravitational waves could be dissipated into the surrounding plasma leading to the 
observed fireball af a GRB. This might be an alternative to the explanation that these 
fireballs are fueled purely by the neutrino flux from the merger, which has the compli- 
cation that the continuous neutrino flux is already much larger than can be explained 
by the present models. 

More work on this subject has to be done, though, before one can make any bold 
statements about the true importance of plasma waves generated by gravitational waves 
in the surroundings of merging neutron stars, leading to gamma ray burst fireballs. Ob- 
viously, the calculations should be extended to spherical symmetries, dipolar magnetic 
fields and spherically decaying gravitational waves. The emws calculated this way will 
probably look more like Bessel functions than plane waves. Still, the results of the first 
order calculations presented in this thesis seem significant enough to motivate further 
research. 
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In this thesis, several theoretical coupling effects between gravitational waves and elec- 
tromagnetic waves were investigated. Some of these effects will be more significant 
than others in an astrophysical context.]] The thesis started with an estimate of the 
general coupling between EMWs and GWs in an EM background field, through the Ein- 
stein field equations. The result of this, and of the subsequent exact calculation, was 
that the coupling efficiency depends on the square of the background field strength and 
the size of the interaction region, with however, and extremely small coupling constant 
(~ 10" 50 (Ns 2 /kg m) 2 ). These conversions therefore occur either in small regions with 
very strong em background fields or in weaker fields extending over very large distances. 

The astrophysical relevance of this process lies in three areas: it could provide indi- 
rect means to detect gravitational waves, it offers a possible explanation for the small 
fluctuations in the cosmic background radiation and it might prove to fuel the fireballs 
produced by gamma ray bursts. These results are discussed here in some detail. 

8.1 Radio waves from binary mergers and magnetars 

In the vicinity of merging neutron star binaries or quaking supernova remnant neu- 
tron stars surrounded by a vacuum, emws with maximum amplitudes of 10 12 V/m and 
50MV/m respectively, could, in principle, be excited. The frequencies of these waves 
are ~ 5 kHz and ~ 10 kHz respectively, viz long radio waves, the latter of which might 
just be detectable with a space based radio array. 

In the more realistic situation where the neutron stars are surrounded by a thin 
plasma, the light generated by the GWs obeys a slightly different dispersion relation, 
with kf = lo 1 — (Aw) 2 instead of the vacuum plane wave dispersion relation k = lo, 
where lo is the driving frequency of the GW. In other words, emws are still produced by 
the GWs, but there is some dispersion of these waves caused by the presence of a plasma. 
This dispersion shifts the wavenumber of the emws by a small amount and also effects 
the polarization by slightly changing the linear polarization relation, B x = —E y . 

Both of these effects vanish, though, in the strong magnetic background fields under 
consideration. These fields suppress the electron mobility in the plasma and thereby 
also the dispersion caused by these electrons, resulting in (Acj) 2 j 0. What is left are, 
again, plane polarized radio waves with the same amplitudes as before in a vacuum. 

The presence of a plasma does determine, however, whether the generated radio 
waves are able to travel over astronomical distances to the earth. The frequency of, 
in particular the emws coming from the binary mergers, is of the same order as the 
the interstellar plasma frequency. The radio waves will therefore be absorbed by the 
plasma, unless close to the merger, non-linear plasma effects result in socalled photon 
acceleration. This effect might lead to higher harmonics of the EMWs with the same en- 
ergy, that might then be able to overcome the interstellar plasma. Light from magnetars 
could have high enough frequencies to overcome the plasma damping without such an 
additional effect. 

7 A11 of the mentioned couplings are only relevant in astrophysics. 
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From the most optimistic point of view, merging neutron star binaries and supernova 
remnants in our local galaxy and the Virgo cluster might, indirectly, be observable in 
radio waves generated by the large amounts of gravitational energy released in these 
processes. On the other hand, the highly idealized calculations (incident plane GWs 
etc.) have probably led to exaggerated results. More detailed (spherical) calculations 
are needed to examine what remains in more realistic situations. 

8.2 Fluctuations in cosmic background radiation 

Another interesting phenomenon that might be explained by the conversion of GWs to 
emws is the fluctuations in the 2.7 K microwave background radiation (grb). 

This radiation results from the decoupling era when the universe was approximately 
20 x 10 4 years old. The temperature had then dropped to 2700 K resulting in the 
recombination of the until then, ionized matter. Because of this recombination, the 
diffusing 'electron mist' evaporated and became transparent to photons, leading to a 
present day, constant flux of thermal photons from this event. 

The wavelength of these photons has reddened because of the expansion of the uni- 
verse during their travel to earth and from this red-shift, one can determine that the size 
of the universe has increased with a factor 1000 since the decoupling. As the present 
size of the universe is supposed to be 10 27 cm (from the redshift of quasars), one finds 
that it was ~ 10 24 cm when the CBR was created. 

If, at that time, a primordial magnetic field of only 10~ 2 Gauss prevailed throughout 
space, a fraction of up to 10 -3 of the GW energy could have been converted into emw 
energy. According to popular cosmological models, this would be enough to explain the 
observed relative fluctuations of 10~ 5 in the cosmic background radiation, which would 
obviously be a very interesting result. 

8.3 Gamma ray bursts 

The most promising candidate as a source for socalled gamma ray bursts (grbs) are 
merging neutron star binaries. These events can produce the required energy to fuel 
GRBs and moreover, they satisfy the temporal and size constraints set by the observa- 
tions. Also, these mergers can be observed as often as once per day, which is comparable 
to the number of detected GRBs. 

A problem in most GRB models is that the energy released by merging binaries is 
released primarily in GWs and not in em radiation. The conversion of the energy in GWs 
to em energy would therefore be very interesting. As was discussed extensively in this 
thesis, the direct conversion of GWs to light is not very effective, and more importantly, 
results in long radio waves, certainly not in gamma radiation. 

A promising alternative is offered by the excitation of magnetohydrodynamic plasma 
waves, such as slow or fast magneto- acoustic waves and Alfven waves. In Section |7.4| , it 
was derived for the first time that, in particular, fast magnetosonic waves can indeed be 
generated through the interaction with gravitational waves passing through the plasma. 
The effectiveness of this process is proportional to the square of the strong magnetic 
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background field and the radio frequency of the perturbing GWs. In Appendix ^| it 
is shown, following |1C|, that nondinear GW effects can cause longitudinal Alfvendike 
waves. 

The magnetosonic waves are generated only when there is a wavelength (or, in the 
static case, a frequency) mismatch between the GW and the excited plasma waves. When 
this is the case, energy from the GW can dissipate into the plasma, which can later on 
emmit the energy in the form of EM (maybe gamma-) radiation. 

Again, the work done in this thesis is not sufficient to fully describe what is going 
on in the vicinity of GRBs, but the result that a fraction of the GW energy from these 
sources could be converted into the observed EM radiation is one that deserves further 
research. Such research would, as a first improvement, have to include dipolar magnetic 
fields and spherical GW solutions for the merging neutron stars. 
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A Constants for exact GW to EMW calculation 

Region I 

Ax = (125) 
+ \ - Bp^fn + e 22 ) - 4°>(&o - 62)} 

d = 

= ( e - 2ifeL -l)[ 

+ 1 { ^ + flf (di + e 33 ) - £l 0) (£00 - 63) } 

+ ^{4 0) 63 + 4 0) ^0-4 0) 6o}] 

Region II 

Ajj = \{^+Bf\i ll+ i 22 )-Ef\i m -i 22 )) (126) 

+ H t + ^ o>6s + 40)60 - 40)60 + 40)62 - 40)63 

Bn = -\{^-Bf\^ ll+ i 22 )-E^{i m -i 22 )) 

- + 4 0) e 30 - sptio + 4 0) 6 2 + b^ 23 
cn = J{^-4 0) (^+^3)-^ 0) (eoo-e33)} 

+ i{§ +^13+^10-^20} 

= -J{^ + 4 0) (^ + ^3)-^ 0) (eoo-e3 3 ) 

- ^{4 0) ei3 + 4°)eio-4 0) 6o} 



Region III 



A "< = fk (127) 
Bin = 

Dm = 
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B Symmetrizing the Maxwell equations 

To prove that {e 5a T P Sa + e^ s W a Sa )B 1 = e^T a s ^B a write: 

e Sa ^T P Sa B^ + e^T a Sa B y - e^T a S7 B a = (128) 

( v aa e^ 5 + v ^e a ^ s -rP^ a8 )Y ff8a B 1 = 

So one has to prove that the expression in brackets is either symmetric in a5 (since T a $ a 
is skewsymmetric in these indices) or vanishes. 
Distinguish between two cases: 

Case 1: For a / ^ 7 the first term gives a = 5 and the other terms vanish, 

Case 2: For a = 7 7^ the second term is zero and the other terms cancel. 

Conclusion : 

( e 8on V P Soi + ^*rr« fa )J3 7 = e^T%B a (129) 

which is nice, because it makes the inhomogeneous Maxwell equations symmetric with 
respect to the homogeneous ones. 
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C Non-linear effects 

It is to be expected that close to a source that emits strong gw's such as magnetars 
and neutron star binaries, non linear GW effects might become important. Consider, 



following 1C], a plane fronted parallel, linearly polarized GW: 



ds z = -dt 2 + a(z - t) 2 dx 2 + b(z - t) 2 dy 2 + dz 2 (130) 

with z — t = u and ab uu + a uu b = 0. Using the same procedure as before, introduce 
an onf, which in this case is a canonical Lorentz frame: 



e (0 j = (1,0,0,0)) (131) 
e (1) * = (Co" 1 , 0,0) 
e {2) * = (0,0, b-\0) 
e (3) * = (0,0,0,1) 

In this metric the gravity induced charge densities do not vanish, and the induced 
current densities have longitudinal components: 



p E = -(lnab) u E 3 (132) 

p B = - (In ab) u B 3 

j E = -(InbUE 1 - B 2 )ex - (lna) u (E 2 + B x )e 2 - (lnab) u E 3 e 3 

3b = -(lnb) u (E 2 + J B 1 )e 1 + (In a) u {E l - B 2 )e 2 - (In ab) u B 3 e 3 

+ V-(nv) = -n(\nab) u (l - v 3 ) 



dn 
~dt 



4- + v-V)v = — (E + vxB) (133) 
at I m 



+ ((lno) u «iei + (ln6) u v 2 e 2 )(l - v 3 ) + ((\na) u v 1 + (lnb) u v 2 )e 3 

If we assume that the GWs are approximately periodic and still have small ampli- 
tudes, the logarithmic factors that appear in these equations can be approximated by 
a ( u ) = S^-oo ex P ^ nam an d b{u) = Y^=-oo bn expinum. To second order in h we 
find (see |Mj for details): 

(lno)„ = iuhe^-V - ^iuh 2 e 2iu1 ^ + c.c. (134) 
(lnb) u = -iuhe^-V - Uuh 2 e 2iw ^ + c.c. 



(lna&) u = — ioj 
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For a GW moving to an initially unperturbed plasma with Eq = Bq = dno/dt = 0, 



(76) and (|133|) , to first order, reduce to: 

dE z 

at 



(In ab) u E z - 47rqn v z (135) 



dv z q 

= —E z and 
ot m 

= §- t (lnab) u E 

Obviously, the flat space plasma oscillations are changed by the GW, especially at 
the resonance frequency oj p = uj, where the parametric excitation of longitudinal plasma 
oscillations is given by: 

^ = -Uujh 2 ^ z E* (136) 
dt 2 y J 

with a growth rate of: 

T = ^uj\h 2 \ (137) 
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